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ABSTRACT 

Soft gamma-ray repeaters (SGRs) and anomalous x-ray pulsars (AXPs) 
are young and radio-quiet x-ray pulsars which have been rapidly spun-down to 
slow spin periods clustered in the range 5 — 12 s. Most of these unusual pulsars 
also appear to be associated with supernova shell remnants (SNRs) with typical 
ages < 30 kyr. By examining the sizes of these remnants versus their ages, we 
demonstrate that the interstellar media which surrounded the SGR and AXP 
progenitors and their SNRs were unusually dense compared to the environments 
around most young radio pulsars and SNRs. We explore the implications of this 
evidence on magnetar and propeller-based models for the rapid spin-down of 
SGRs and AXPs. We find that evidence of dense environments is not consistent 
with the magnetar model unless a causal link can be shown between the 
development of magnetars and the external ISM. Propeller-driven spin-down by 
fossil accretion disks for SGRs and AXPs appears to be consistent with dense 
environments since the environment can facilitate the formation of such a disk. 
This may occur in two ways: 1) formation of a "pushback" disks from the 
innermost ejecta pushed back by prompt reverse shocks from supernova remnant 
interactions with massive progenitor wind material stalled in dense surrounding 
gas, or 2) acquisition of disks by a high velocity neutron stars, which may be 
able to capture a sufficient amounts of co-moving outflowing ejecta slowed by 
the prompt reverse shocks in dense environments. 

Subject headings: Stars: neutron — pulsars: individual (SGRs, AXPs) — ISM: 
supernova remnants 
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1. Introduction 



Soft gamma-ray repeaters (SGRs) are neutron stars whose multiple bursts of 
gamma-rays distinguish them from other gamma-ray burst sources (e.g [Hurley 2000| for a 
recent review). SGRs are also unusual x-ray pulsars in that they have spin periods clustered 
in the interval 5 — 8 s, and they all appear to be associated flCline et al. 1980] ; |Felteu 1981 



Kulkarui fc Frail 1993j ; |Vasisht et al. 1994| ; |Hurley et al. 19994 |Hurley et al. 1999cj ; |Corbel 



et al. 1999 ; Cliue et al. 2000 ) with supernova remnants (SNRs), which limits their average 
age to approximately 20 kyr ( |Braun, Goss, fc Lyne 198^ ). The angular offsets of the SGRs 
from the apparent centers of their associated supernova remnant shells indicate that SGRs 
are endowed with space velocities > 500 km s _1 , which are greater than the space velocities 
of most radio pulsars flCordes fc Chernoff 1998]) . Anomalous x-ray pulsars (AXPs) are 
similar to SGRs in that they are radio quiet x-ray pulsars with spin periods clustered in 
the range 6 — 12 s, and have similar persistent x-ray luminosities as the SGRs (~ 10 35 
ergs s -1 , see e.g. [Stella, Israel, fc Mereghetti 1998| for a recent review). Most of the AXPs 



appear to be associated with supernova remnants, and therefore they are also thought 
to be young neutron stars like the SGRs. The spin periods of both AXPs and SGRs are 
increasing with time (spinning-down), and show no evidence for intervals of decreasing spin 
period (spin-up), although the spin-down rates of many of the SGRs and AXPs appear to 
be variable or "bumpy" (e.g. Baykal fc Swank 1996 ; Marsden, Rothschild, fc bingcnlcltcr 
[Woods et al. 1999c| ; |Woods et al. 2000| ; although for a different viewpoint see |Kaspi,| 



1999 



Uhakrabarty, fc Steinberger 1999|). 



The lack of identified companions at non x-ray wavelengths (e.g. [Stella, Israel, fc 



Mereghetti 1998|) and Doppler shifts associated with binary orbital motion QMereghett 



Israel, fc Stella 1998| ), together with the problem of accelerating binaries to space velocities 
> 1000 km s" 1 ( [Rothschild, Kulkarni, fc Lingenfelter 1994Q , imply that SGRs and AXPs 
are not members of high mass binary systems, although low mass systems with neutron star 
companions of < 1M are not constrained in most cases. If SGRs and AXPs spin-down 
primarily via the emission of magnetic dipole radiation (MDR), as do radio pulsars, 
then they must have surface dipole fields of ~ 10 14 G or greater (i.e. "magnetars"; e.g. 



Thompson fc Duncan 199~5| ). Observations of the SGRs 1806-20 (|Kouveliotou et al. 1998|) 
and 1900+14 ( |Marsden, Rothschild, fc Lingenfelter 1999| ; |Woods et al. 1999~c|) , however, 
indicate that the present-day spin-down rates of these SGRs are inconsistent with simple 
MDR, given the ages of their associated supernova remnants ( Harding, (Jontopoulos, fc] 
Kazanas 1999| ; Rothschild, Marsden, fc Lingenfelter 1999]) , and imply that the spin-downs 
are due to winds. Magnetar-strength fields might still be possible in these sources, however 
if Alfven wave wind emission is infrequent and intermittent ( |Harding, Gontopoulos, fc 



Kazanas 1999| , so that the presently observed spin-down rates are atypical. Alternative 
scenarios for SGRs and AXPs involving typical pulsar magnetic fields (~ 10 12 G) have been 
proposed ( |van Paradijs et al. 1995] ; Alpar 2000; Chatterjee, Hernquist & Narayan 2000; 
Chatterjee & Hernquist 2000). In these models, the SGRs and AXPs have spun-down 
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rapidly via magnetospheric accretion torques from outflowing "propeller effect" winds. The 
assumed sources of the magnetospheric material are either fallback accretion disks (Alpar 
2000; Chatterjee, Hernquist & Narayan 2000; Chatterjee & Hernquist 2000), or fossil disks 



formed from expanding supernova ejecta intercepted by high velocity neutron stars ([van 
Faradijs et al. 1995| ; |Uorbet et al. 1995|) . 



Here we present a fresh look at evidence which suggests that SGRs and AXPs are 
born into unusually dense environments. We show that the environments of the SGR and 
AXP progenitors into which their SNRs expand are the dense, warm and cool phases of 
the interstellar medium (ISM), and not the hot tenuous phase of the ISM where most of 
the neutron-star-producing, core collapse supernovae of massive O and B stars occur and 
where most young radio pulsars are found. This implies that there is an environmental 
factor influencing the development of SGRs and AXPs. The structure of this paper is 
as follows. We first discuss the typical environments of supernova progenitors in §2, and 
then supernova remnants associated with SGRs and AXPs in §3, followed by a discussion 
of the SGR and AXP ages and distances are in §4. In §5 the density of the SGR and 
AXP progenitor environments is discussed, and in §6 a similar analysis is done for the 
SNRs associated with young radio pulsars. In §7 we discuss the statistical significance of 
the results and the implications for magnetar and propeller-based models for SGRs and 
AXPs. Finally, the Appendix contains a short discussion of the ages, distances, and other 
information for each SGR and AXP. 



2. The Environments of Supernova Progenitors 



What are the environments typical of neutron star progenitors? Observations clearly 
show that the majority of neutron stars are formed in "superbubbles" : evacuated regions 
of the ISM which surround the OB associations in which the massive progenitors of most 
neutron stars live. This is because most O and B stars (> 80%) are observed (e.g. McCray 
& Snow 1979) to occur in clusters formed from giant molecular clouds (> 1O 5 M ); much 
smaller clouds are disrupted by the radiation and winds from the first O star that forms. 
These massive (> 8 M : |Woosley fc Weaver 1995[ ) and slow-moving (~ 4 km s _1 : |Blaauw 



1961| ) O and B star progenitors of Type II and Ib/c supernovae do not travel far from 



their birthplaces during their relatively short (< 30 Myr: (Schaller et al. 1992 ) lives. The 



supernovae from these massive stars are therefore heavily clustered in space and time and 
form vast (> 100 pc) HII regions/superbubbles (e,g. [Mac Low fc McCray 19881) filled with 
a hot (> 10 6 K) and tenuous (n ~ 10~ 3 cm -3 ) gas. 

There is also direct observational evidence that the great majority (>80%) of neutron 
star stars are born into superbubbles consisting of hot and diffuse ISM gas. Observations 
of a sample of 49 spectroscopically identified Type II and Ib/c, core-collapse supernovae 
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occurring in face-on late-type spiral galaxies by van Dyk et al. (1996) found that 72% ±10% 
of Type II and 68% ± 12% of Type Ib/c supernovae are in resolvable giant HII region 
superbubbles. We suggest that these fractions are, in fact, only a lower limit on the 
occurrence of core-collapse supernovae in superbubbles, because of the difficulty in detecting 
faint HII regions in distant galaxies. Correcting for this effect using the Ha luminosity 
distributions of Kennicutt, Edgar, & Hodge (1989), we find that the Ha threshold used by 
van Dyk et al. (1996) would have allowed them to resolve only 76% of the HII regions in 
their sample galaxies. This clearly suggests that the great majority - i.e. 90% ± 10% - of 
neutron star progenitors reside in the hot and diffuse {n < 0.01 cm -3 ) ISM in superbubbles, 
with < 20% occurring in the denser phases of the ISM. 

When one corrects for selection effects, the distribution of surface brightness of 
all detected Galactic supernova remnants provides another measure of the fractions of 
supernovae occurring in the warm and hot phases of the ISM. Detections of radio remnants 
are limited by surface brightness, which is independent of their distance, and SNRs in the 
denser ISM have larger surface brightnesses than SNRs of the same age in the hot / diffuse 
ISM (e.g. |Gull 19731) . A study flHigdon fc Lingenfelter 19~80| ) of the age versus surface 



brightness of the remnants of historical supernovae, using Gull's model, suggests that the 
maximum detectable ages of radio supernova remnants above a nominal surface brightness 
density of 10~ 21 W m -2 Hz -1 sr _1 at 408 MHz are roughly 20 kyr in the warm (n = 0.1 
cm~ 3 ) gas and 6 kyr in the hot superbubble (n = 0.001 cm -3 ) gas. Combining this with 
an analysis of the number distribution of SNRs as a function of surface brightness suggests 
that about 30% of the detected radio SNRs are in the warm denser (n > 0.1 cm -3 ) phase 
of the ISM (e.g. [Higdon fc Lingenfelter 1980|) . This implies that only about 10% of the 



supernovae occur in the warm (n > 0.1 cm -3 ) ISM 0, and about 90% in the hot phase, 
since, for a total (including Type la) SN rate of 1 SN every 30 yrs, there are then roughly 
180 SNR younger than 6 kyr that are detectable in the hot superbubble ISM, and roughly 
70 SNR younger than 20 kyr that are detectable in the denser warm ISM, or about 30% of 
the detected radio SNRs in the warm ISM - as the observed distribution implies. 



3. The Supernova Remnants Associated With SGRs and AXPs 



The environments of SGRs and AXPs are probed by the blastwaves of their associated 
supernova remnants, and from the size of the remnant shell as a function of the age we can 



2 We should mention that Higdon & Lingenfelter (1980) further suggested that this 
fraction represented the filling factor of the warm ISM, on the mistaken assumption that 
the supernovae were uniformly distributed, rather than being predominantly clumped in OB 
associations. 
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constrain the external density. In Table 1 we have listed the 12 known SGRs and AXPs 
and their associated supernova remnant shells. For more information on the individual 
SGRs and AXPs, see the recent reviews of Hurley (2000a) and Mereghetti (1999), and the 
Appendix. We include the source AXP 1845-0258 - although no period derivative has been 
measured for this source - because it is commonly cited as an AXP (see e.g. the discussion 
Gaensler et al. 1999| ). We also include the AXP 0720-3125 for completeness, although 



in 



its identification as an AXP has been questioned as well (see Appendix). The new and 
tentative SGR 1801-23 has also been included on the basis of two soft gamma-ray bursts 
consistent with a single location on the sky (|Gline et al. 



The identification of the associated remnants are based on positional coincidences 
between the remnant and the SGR/AXP, and in some cases on similar distances as implied 
by the hydrogen column densities measured from the x-ray spectrum of the SGR/AXP and 
its associated remnant. Based on these criteria, associated SNRs have been previously been 
identified for 7 out of the 12 SGRs and AXPs (see Table 1 and Appendix). Using these 
same criteria, we suggest three additional SNR associations with SGR/AXPs (also included 
in Table 1 and Appendix). We suggest two new AXP/SNR associations, AXP 1709-40 
with the SNR G346.6-0.2 QGreen 2000| and references therein) and AXP 1048-5937 with 



G287.8-0.5 ( |Jones 1973|) , based on the near coincident positions between each AXP and a 



known supernova remnant, and also (for AXP 1048-5937/G287.8-0.5) on similar implied 
distances for the AXP and the remnant. These two associations may have been discounted 
previously because they would imply larger than average neutron star velocities, but in 
view of the comparable velocities implied for three of the SGRs these possible associations 
should be considered. AXP 1048-5937 lies 2.2 shell radii from the apparent center of the 
SNR G287.8-0.5, which is associated with the Carina Nebula, a region of massive star 
formation at a distance ~ 2.5 — 2.8 pc ( The fc Vlccming 1971 ; Seward fc Chlcbowski 



1982| ). This relatively small distance of the SNR is quite consistent (see Fig. 2) with the 
low Nh = (0.45 ± 0.10) x 10 22 cm" 2 found from the AXP 1048-5937 x-ray spectral fits 
( poster broeck et al. 1995| ). For an estimated SNR age of ~20 kyr, the displacement of 
the AXP implies a relatively modest neutron star velocity of ~ 1000 km s _1 . A similar 
transverse velocity is indicated by the association of AXP 1709-40, which lies 1.7 shell 
radii from the center of the well defined shell remnant G346.6-0.2, assuming a similar 
age and a distance of 3 to 5 kpc, based on Galactic structure. Such a SNR distance is 
quite consistent with its association with AXP 1709-40, which has an x-ray absorption 
N H = (1.81 ± 0.07) x 10 22 cm" 2 ( gugizaki et al. 19971) . We also propose that the new 
SGR candidate 1801-23 ( |Gline et al. 2000| ) may be associated with the well studied SNR 



W28. Although there is no distance estimate to the SGR flCline et al. 2000| ), W28 is the 
only known SNR through which the thin SGR error box passes. Therefore we regard the 
tentative W28/SGR 1801-23 association as promising and encourage observations to search 



for x-ray point sources (e.g. [Andrews et al. 1983j ) within the remnant which may be the 



SGR counterpart. We include these three tentative SNR associations because they are 
suggestive, but as we show in §7.1 these data can be left out of the sample without affecting 
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the overall conclusions of our analysis. 

Using the SNR catalog of Green (2000), we are unable to identify any possible remnants 
associated with AXPs 0720-3125 and 0142+615, which may be the two oldest AXPs based 
on their pulsar timing ages 0.5P/P, where P is the pulsar period and P is the period 
derivative, of 40 kyr ( [Haberl et al. 1997] ) and 60 kyr (e.g. [White et al. 1996|) , respectively. 
For AXP 0720-3125 the lack of an associated remnant is not surprising, because this source 
is thought to be so close (~ 100 pc) that a neutron star with a velocity of ~ 1000 km s _1 
would have traveled ~ 40 pc in 40 kyr and hence could have originated from a large area 
of the sky. In addition, the large scale radio surveys of the Galactic Plane appear to be 
incomplete for very large (> 1°) diameter remnants (e.g. Duncan et al. 19"9~7| ), indicating 



that such a remnant could easily go undetected. AXP 0142+615 is situated in or behind a 
molecular cloud complex ( [Israel, Mereghetti, fc Stella 1994] |White et al. 1996]) , and if its 



associated remnant expanded in high density material it may have already passed into the 
radiative phase and faded below the surface brightness detection threshold. Therefore the 
detection of the remnants associated with both AXP 0720-3125 and AXP 0142+615 would 
be difficult, and we can not assign meaningful limits to the physical size of their (unknown) 
associated remnants, given the present data. We encourage new deep observations of the 
regions surrounding these objects to look for associated supernova remnants. 

We see from Table 1 that 10 out of the 12 SGR/AXPs have proposed associations 
with radio shell SNRs. The probabilities of a chance coincidence for some of the individual 
associations have been estimated (e.g. |Cline et al. 1980] ; [Felten 1981| ; |Kulkarni fc Frail 



1993|) at a few times 10 2 or less (see Appendix). The chance probability for an association 



between an SGR/AXP and an SNR can be estimated by simply considering the spatial 
density of known supernova remnants in the Galactic plane (in the following we assume 
that the localization error box of the SGR/AXP is small compared to the size of the SNR). 
For the extreme case, we consider just the SGRs and AXPs located in the highly crowded 
inner Galaxy. There all 8 SGR/AXPs appear to have associated SNRs. These have Galactic 
longitudes between 286° and 43° and latitudes within ± 1.2°. Within this area of 1.0 xlO 6 
arc min 2 , there are 142 known SNRs (Green 2000) covering a combined surface area of 
3.4xl0 4 arc min 2 , which is actually an overestimation due to overlap between the remnants. 
For a single association with the SGR or AXP displaced from the center of the remnant 
by as much as 2.3 times the remnant radius (as implied by the SGR 1627-41/G337. 0-0.1 
association), the chance association probability is roughly 2.3 2 x 3.4 x 10~ 2 ~ 18%, or 
~ 16% correcting for overlap. Thus there is a significant probability that a single association 
might be simply a chance coincidence, if we relied solely on position, but the probability of 
getting even 4 chance associations out of 8 tries is only 2%, and the probability that 8 out 
of 8 are spurious is just 4 x 10~ 7 . 



4. The Distances and Ages of the Associated Remnants 
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In order to probe the SNR environments, we need to know both the radii of the 
associated SNR shells, which can be determined from their measured angular diameters, 
and their estimated distances. These distances are listed in Table 2, and discussed in 
the Appendix. All but three of the distances were determined from Galactic kinematic 
arguments based on interactions or associations with adjacent objects with known distances 
(e.g. HII regions or molecular clouds), or on absorption line measurements. We did not 
use distances estimated from the standard SNR surface brightness/diameter relationships, 
because we did not want to bias the SNR sizes toward SNRs in the dense ISM, which is 
where most of the observed radio SNRs are located (e.g. Higdon & Lingenfelter 1980; 
Kafatos et al. 1980). SNRs expanding in the tenuous hot ISM become large and hard to 
detect (e.g. Puncan et al. 1997 ) much quicker than SNRs in dense environments and are 



under-represented. 

The positions of the Galactic SNRs associated with SGRs and AXPs are shown 
in Figure 1, projected onto a model of the Galactic plane distribution of free electrons 
( [lay! or fc Gordes 1993 ) which traces star formation regions delineating the spiral arms 



(e.g. Georgelin & Georgelin 1976). We see that the estimated distances of the SNRs are 
quite consistent Galactic structure, placing them in or near the spiral arms where the bulk 
of recent massive star formation has occurred. We also see that the estimated distances 
of the SNRs are not systematically underestimated in our analysis, since 5 out of the 9 
SNR-AXP/SGR associations lie roughly at or beyond the distance of the Galactic center - 
as would be expected for a relatively unbiased sample of Galactic sources. This also shows 
that they are not likely to be systematically much farther away, or most would lie on the 
other side of the Galaxy or even outside of the Galaxy. 

To further test both the AXP/SGR-SNR associations and their estimated distances, 
we compare the electron column depth associated with the SNR distance with the H column 
depth determined from the x-ray absorption of the AXP/SGR emission. We calculate the 
dispersion measure to each source given the assumed distances and compare that to the 
measured Nh value for each source. Table 2 lists the distance ranges, calculated dispersion 
measures using the Taylor & Cordes (1993) model, and observed Nh values for each 
SGR/AXP and associated SNR. The dispersion measure (DM) vs. Nh is plotted in Figure 
2. As seen from the Figure, the data points are well correlated between lines representing 
free electron to hydrogen ratios e/H = 0.03 — 0.07, which is consistent with the range of 
0.02 — 0.08, expected from the mean mass fraction of ionized gas in the Galaxy of ~ 1 — 5% 
(e.g. Spitzer 1978] ). Thus we see that each of associations is consistent both in proximity of 



direction and in similarity of distance. This also shows that the estimated distances to the 
SGRs/AXPs and their associated remnants are not systematically underestimated. More 
information on the Nh values for each source is given in the Appendix. 

The ages of the most of the SGRs and AXPs are much more uncertain than their 
distances. This is because the usual pulsar timing age formula - appropriate for pulsar 
spin-down due to magnetic dipole radiation (MDR) torque (although |Gaensler fc .brail] 
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2000) - does not seem to be applicable for SGRs and AXPs, in general, because it yields age 
estimates which are inconsistent with the likely ages of the associated supernova remnants 
and the implied neutron star velocities for at least three of the sources (SGR 1806-20, 
SGR 1900+14, and AXP 2259+58.6: [Kouveliotou et al. 199Sj ; |Marsden, Rothschild, fc 



Lingenfelter 1999| and |Gorbet et al. 1995| , respectively). In addition, all of the SGR and 
AXPs for which pulsar braking indices have been measured - AXPs 1709-40 and 2259+58.6 
( |Kaspi, Ghakrabarty, fc Steinberger 19^9] ) and SGR 1806-20 QWoods et al. 2000| ) - have 
braking indices which are inconsistent with simple MDR spin-down. The MDR timing age 
for AXP 1841-045 of 4000 yr, however, is consistent with the estimated age of the associated 
SNR Kes 73 flGotthelf, Vasisht, fc Dotani 1999|) , but the braking index for this AXP has 
not yet been measured. We therefore adopt ages for the SGRs and AXPs based on the 
ages of their associated SNRs. Some of the associated remnants (N49, Kes 73, G29. 6+0.1, 
and CTB 109) are relatively well-studied and have age estimates based on observed shock 
velocities and/or x-ray temperatures. For these remnants we take the minimum and 
maximum values for the ages of each remnant from the literature. For the rest of associated 
SNRs, we adopt a lower age limit of U ow = min[tf e ,t v ], where tf e = d m i n 9sNR/v e j is the 
age of the remnant assuming free expansion, with d m i n the minimum estimated distance, 
6snr the SNR radius in radians, and v e j ~ 10 4 km s _1 is the maximum ejecta speed in free 
expansion. In addition, if we assume that the transverse velocity of the associated SGR or 
AXP cannot exceed v max , the minimum age is t v = d min 9*/v max , where 9* is angular offset 
of the pulsar from the center of the remnant. In this paper we will assume v max = 2000 km 
s _1 , which exceeds the observed velocity of any pulsar (e.g. |Cordes fc Chernoff 1998|) . For 
the upper age limit of the SNRs we choose 30 kyr, which is the maximum estimated age for 
the SNR/radio pulsar associations shown in Table 30. For each SNR, we assume that the 
most likely age is the arithmetic mean of the age range, which is equivalent to assuming 
that the supernova rate is uniform in time throughout the Galaxy. 



5. The Densities of SGR/ AXP Progenitor Environments 



In Figure 3 we have plotted the SNR shell radii Rsnr as a function of the estimated 
age t of each remnant associated with an SGR or AXP. Overplotted with solid lines are 
simple approximations of the evolutionary tracks ( |5hull, Fesen, fc Saken 1989|) of supernova 
remnant expansion in a wide range of the external ISM densities n, assuming a total 
supernova kinetic energy of 10 51 ergs. These SNR evolutionary tracks move through three 
phases: the initial free expansion/ejecta-dominated phase of the remnant, where the mass 



3 The overall conclusions are unaffected if we instead choose 20 kyr ( [Braun, Goss, fc Lyne 
1989|) for the maximum observable SNR age, as discussed in § |7.1.| 
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of the SN ejecta is much greater than the mass of the swept-up ISM and Rsnr oc t; the 
Sedov/adiabatic phase, which begins when the mass of the swept-up ISM is roughly > 10% 
of the mass of the SN ejecta, and the remnant slows down with Rsnr oc t 2 / 5 ; and finally 
the radiative/snowplow phase, where the shell of swept-up ISM radiates away the energy 
of the remnant and it slows further with Rsnr oc t 2 / 7 (e.g. [Shull, Fesen, fc Saken 1989|) . 



For individual supernova remnants, the tracks may differ by less than ±15% for a factor of 
three variation in supernova kinetic energies ( |Woosley fc Weaver 1995| ) . Also overplotted as 



dotted lines are the tracks of neutron stars born at the origin of the supernova explosion 
with varying velocities, showing the times required for them to catch up with the outer 
(radio) shell of the expanding remnant. 

We see from Figure 3, that in spite of the uncertainties in the remnant ages and 
distances, all of the supernova shell remnants associated with SGRs and AXPs seem to 
reside in the denser (n > 0.1 cm -3 ) phases of the ISM, where < 20% of the neutron star 
producing supernovae occur. Although the ages of the the SNRs are highly uncertain, 
information on the progenitor ISM density can be extracted from the SNRs because 
densities inferred from the Sedov and Radiative phase expansion formulae depend much 
more strongly on the SNR radii than on their ages. As we discuss in more detail below, 
the inferred distribution of progenitor ISM densities is not consistent with an origin of 
SGRs and AXPs that results from a purely intrinsic property of the neutron stars, because 
such stars should predominantly be born in the hot, diffuse (n ~ 10~ 3 cm" 3 ) phase of the 
ISM, where > 80% of the neutron star and pulsar producing SN occur. The probability of 
finding 10 such SNRs out of 12 possible occurring in the denser {n > 0.1 cm -3 ) phases of 
the ISM, when at most ~ 20% are expected is only ~ [12!/(10!2!)](0.2) 10 (0.8) 2 ~4x 10" 6 . 
Even if we assume the maximum 16% probability for chance associations with SNR for 
all SGRs and AXPs and assume that all chance associations will be in denser ISM, the 
combined probability of any single, typical SNR being in the denser ISM is 0.33, equalling 
0.2 naturally occurring plus 0.8 x 0.16 occurring by chance. Thus the probability of finding 
10 out of 12 possible occurring in the denser ISM is still only ~ [12!/(10!2!)](0.33) 10 (0.67) 2 
or ~ 4 x 10" 4 . 

There are two obvious ways that the distribution of these SNRs might be made 
consistent with that of typical neutron stars and pulsars, i.e. 80% residing in the hot diffuse 
ISM and 20% in the dense ISM. One way would be if the SNR distances were systematically 
underestimated by roughly a factor of 3, but as we saw in Figure 1 and 2, that is not 
consistent with Galactic structure and would place roughly half of the SNRs well outside 
the Galaxy. Alternatively, exceedingly weak (E << 10 51 ergs) explosive energies for the 
SGR and AXP supernovae could in principle explain the unusually compact remnants 
associated with SGRs and AXPs, but this is inconsistent with dynamical requirements of 
supernova explosions (e.g. Woosley & Weaver 1995), and with the direct observations of 



SGRs and AXPs associated SNRs (e.g. |Vancura et al. 1992| ; |Rho fc Petre 1997| ; |Gotthell 



fc Vasisht 1997| ; [Long et al. 1991| ). Moreover, there is direct evidence that at least 8 out 
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of the 10 SGRs and AXPs associated SNRs are indeed in dense environments. OH maser 
emission, attributed to SNR shock interactions with molecular clouds, has been detected 
(Frail et al. 1996; Green et al. 1997) from molecular clouds thought to be associated with 
five of the SNRs (CTB 33, Kes 73, W28, G10.0-0.3, and G346.6-0.2). CTB 109 (e.g. Huang 
& Thaddeus 1985; Tatematsu et al. 1987), G287.8-0.5 QJones 1973| ), and N49 QHughes et 



al. 1984 ) also show evidence for molecular clouds associations. These associations with 



molecular clouds clearly show that the supernova remnant shells associated with SGRs and 
AXPs are expanding in high density environments. 



6. Comparison of SGRs/ AXPs to Young Radio Pulsars 



We can compare the observed properties of the SNRs associated with SGRs and AXPs 
with the remnants associated with another population of young neutron stars - the young 
radio pulsars listed in the catalogs of Taylor, Manchester, & Lyne (1993) and Taylor et al. 
(1995). We consider only the youngest pulsars with timing ages less than ~ 30 kyr, so as to 
be comparable to the age range in the SGR/AXP sample. The MDR timing ages of radio 
pulsars, unlike those for the SGRs and AXPs, are generally though to be good measures 
of their true ages (e.g. Gordes fc Chcrnoff 1998|) , and we assume that the remnant ages 



are consistent with the MDR timing ages of their associated pulsars. A notable exception 
is PSR J1801-2451, which has an estimated age much greater than its MDR timing age 
( Gaensler fc Frail 2000|) - possibly indicating the presence of non-MDR spin-down torques 
in this pulsar. For the pulsar distances, we use the distances derived ( Taylor, Manchester 



|fc Lyne 1993| ; lay lor et al. 1995 ) from the measured dispersion and the Taylor & Cordes 



(1993) model of the free electron distribution in the Galaxy. 

The supernova remnant shells of many of the young pulsars in Table 3 have not been 
detected in large scale radio surveys of the Galactic Plane (see e.g. [Whiteoak fc Green 



199B1 |Reich et al. T98g |D uncan et al. 19971 ), or in deep radio observations of the fields 



surrounding young pulsars (e.g. |Braun, Goss, fc Lyne 1989| ; [Frail, Goss, fc Whiteoak 1994j ). 



For these pulsars, we found no likely associated SNRs in the Green (2000) catalog within a 
search radius corresponding to a transverse velocity of as much as 2000 km s , assuming 
the estimated distances and timing ages for these sources listed in Table 1. Since the Green 
(2000) SNR catalog is more or less complete ( |Whiteoak fc Green 1996] ) down to a limiting 
surface brightness, the missing SNRs for these objects have probably expanded such that 
their surface brightnesses have faded below the limiting surface brightness of the radio 
surveys of the Galactic plane. We therefore assume that the radio shells of the remnants 
associated with these pulsars (except for the Crab, for which we use a lower limit of 17 
pc for the remnant radius; [Frail, Goss, fc Whiteoak 1994|) have expanded beyond this 
detectability threshold, and following Braun, Goss, & Lyne (1989) we assign lower limits of 
30 pc to the undetected remnants corresponding to these pulsars. The observed and derived 
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parameters for the young radio pulsars are listed in Table 3, with references to the timing 
data, distances, and SNR shell radii. A plot of the supernova remnant shell radii versus age 
for the radio pulsars is shown in Figure 4. 

We clearly see from Figure 4 that most of the young pulsars appear to have been born 
in the hot, diffuse phase of the ISM, as expected from other observations of their O and B 
star progenitors and their environments discussed in the preceding section. Only ~ 5 of 
the 16 young pulsars have SNRs which may be expanding in the denser phases (n > 0.1 
cm -3 ) of the ISM, which is consistent with the expected fraction of ~ 20% or less despite 
the small number statistics. These results are summarized in Table 4. 



7. Discussion 



7.1. The Significance of the Environmental Evidence 



As seen from Figures 3 and 4, the SGRs and AXPs appear to form in denser regions 
of the interstellar medium than the sample of young radio pulsars. The significance of the 
apparent disparity between the SGR/AXP and radio pulsar environments can be evaluated 
by using "survival analysis" ( [Miller 1981| ) methods, which are statistical techniques 
incorporating "censored" data (e.g. upper limits) into data analysis. We first convert the 
age and SNR radii values of Tables 1 and 3 to ambient ISM densities n using the standard 
formulae in Shull (1983) and using values for each age and SNR radius midway between 
the range of values listed in the tables. A total supernova kinetic energy of 10 51 ergs was 
assumed for both the SGR/AXPs and radio pulsars. The lower limits on the radii of the 
undiscovered SNRs in Table 3 therefore became upper limits on the ambient densities for 
these SNRs. 

The resulting distributions of n for the SGR/AXPs and radio pulsars were then tested 
for consistency with a single parent sample using the statistical analysis package ASURV 
Rev 1.2 ([Lavalley, Isobe, fc Feigelson 1992]) , which implements the methods presented 
in Feigelson & Nelson (1985). The two-sample univariate nonparametric tests used to 
compare the distributions consisted of two versions of the Gehan's Generalized Wilcoxon 
test (permutation and hypergeometric variances), the Logrank test, and the Peto & Peto 
and Peto & Prentice Generalized Wilcoxon tests (see [Lavalley, Isobe, fc Feigelson 1992| and 



references therein; [Miller 1981| ; [Feigelson fc Nelson 1985Q . These methods differ in their 



assumptions regarding the censoring process and in the ways which they weight the data, 
but we found that they gave results which were consistent to within two orders of magnitude. 
We use the statistical tests to compute the probability P sam e that the distribution of n for 
the SGRs and AXPs is consistent with the same distribution for the young radio pulsars. 
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For the full sample of SNR associations in Table 1 we calculate P sarne = 9 x 10 5 — 6 x 10 4 . 
If we exclude from the sample the objects with tentative SNR associations - SGR 1801-23 
and AXPs 1709-40 and 1048-5937 - we find P same = 1 x 10~ 4 - 8 x 10" 4 . Finally, if we 
change the maximum possible remnant age from 30 kyr to 20 kyr ( Praun, Goss, fc Lyne 



1989| ), we obtain P sam e = 5x 10 4 — 3 x 10 3 for the full set of SNR associations in Table 



1, and Psame = 8 x 10 — 5 x 10 if we exclude the tentative SNR associations. The 
data therefore supports the conclusion that the progenitors of SGRs and AXPs exploded in 
denser environments than the progenitors of radio pulsars. 

This conclusion can be alleviated only if one assumes that a large fraction of the 
previously claimed SNR associations with AXPs and SGRs are spurious and the remnants 
truly associated with these objects have large diameters and are currently undetected. As 
discussed in §3, however, it is unlikely that more than a couple of the associations in Table 
1 are spurious, given the distribution of supernova remnants in the Galactic plane. At least 
five of the SNR associations in Table 1 would have to be spurious before one could conclude 
that the evidence for dense SGR and AXP progenitor was insignificant (i.e. the probability 
Psame > 0.01). The association of SGR 1900+14 with SNR G42.8+0.6 has recently been 
questioned after the recent discovery ( [Lorimer fc Xilouris 2000]) of a 40 kyr old pulsar 
PSR J1907+0918 equally close ~ 20' to the SNR, and the discovery (|Vrba et al. 2000|) of 
a compact cluster of massive stars at an estimated distance ranging from 12 to 15 kpc, 
and only ~ 10" from the SGR 1900+14 line of sight. The relatively low H column depth 
from absorption in the SGR's x-ray spectrum ( |Hurley et al. 1999 b| ; |Woods et al. 1999"a| ) 
clearly implies (e.g. Figure 2) a distance closer to 5 kpc, so that the SGR and and SNR 
may not be related to the cluster. If SGR 1900+14 is associated with the compact cluster 
at the more extreme distance, Vrba et al.(2000) suggest an association with a possible 
compact (< 1 pc diameter) SNR indicated by the local diffuse x-ray emission. This new 
SGR/SNR association would therefore imply a much greater SGR 1900+14 progenitor 
density than the G42.8+0.6 association, increasing the evidence for dense progenitor 
environments for SGRs. A similar compact (< 1.0 pc diameter) cluster of massive stars has 



also been observed flFuchs et al. 19991) to be coincident with SGR 1806-20 and G10.0-0.3. 
In addition, AXP 1048-5937 and its associated SNR G287.8-0.5 are associated with the 
Carina star- forming region ( [The fc Vleeming 1971| ; |Seward fc Chlebowski 1982] ). In these 



cases the associations with the clustered star-formation regions are quite consistent with 
the SGR/SNR associations, since the distances of the star-forming regions and the SNRs 
are comparable. 



7.2. SGRs and AXPs as Magnetars 



The unusual properties of SGRs and AXPs may be due to the fact that they have 
unusually strong magnetic fields, as in the magnetar model ( |Duncan fc Thompson 1992 
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Tliompson h Duncan 1993: Thompson fc Duncan 1995 ; Thompson fc Duncan 1996 ). 



In the context of this model, it has been suggested that the association of SGRs and 
AXPs with dense progenitor environments may be explained by a selection effect in the 
following manner. Magnetars are thought to form from progenitor stars with high angular 
momenta, such that the superstrong (10 14 — 10 15 G) magnetic field can be generated in the 
protoneutron star by dynamo action just after core collapse (Puncan fc Thompson 1992 ). 



Since the total stellar angular momentum J oc MRvr (assuming rigid rotation, where vr is 
the observed rotational velocity), one might expect the largest and most massive stars to 
be the likely progenitors of magnetars. Furthermore, since the main sequence lifetimes of 
stars decrease with increasing mass, the more massive stars would be the first to supernovae 
- therefore exploding into more dense surroundings before the parent molecular cloud has 
been cleared by the successive supernovae of less massive stars. 

The two basic assumptions of this scenario do not appear to be consistent with 
observations. First, the observed stellar rotation velocities (e.g. prilling fc Landoldt 20001) 



show that the angular velocity Vt increases with stellar mass only up to about 2 M and 
then decreases with increasing mass, dropping to only 20% of its maximum value for stellar 
masses > 25M . Since the stellar evolution timescale even in massive stars is much longer 
than the convective timescale (e.g. |Kippenhahn fc Weigert 1991| ), the rotation of the star is 
effectively rigid and one would therefore expect the pre-collapse cores of the most massive 
stars to have less angular momenta than the cores of less massive supernova progenitors. 
This is supported by detailed calculations of the stellar evolution of rotating massive stars 
up to the time of core collapse ([Heger, Langer, fc Woosley 2000| ), which indicate that the 
final angular momentum of the core region is less for a massive star than for a less massive 
star with the same initial rotation rate. This is due to the greatly increased mass losses of 
luminous massive stars, which shed angular momentum in their powerful stellar winds (e.g. 
Heger, Langer, fc Woosley 2000| ). 



Second, the supernovae of massive stars would still not be expected to all explode in the 
denser medium, because star formation in massive (> 10 5 M ) clouds is episodic through 
several generations (e.g. [McKcc fc Williams 1997| ). These successive generations commence 
roughly every ~4 to 5 Myr, over a period of ~20 to 30 Myr, each producing roughly 10 3 O 
and B star supernova progenitors. Each subsequent generation of star formation commences 
shortly after the onset of supernova explosions from the previous generation, at the end of 
their main sequence lifetimes of >3 Myr. Therefore only the most massive stars in the first 
generation of star formation in a new star forming region would be expected to explode in 
denser environments, and the subsequent generations of massive stars should be born in 
the same low density superbubble environment as the population as a whole. Thus even if 
magnetars were produced by the most massive stars, we still would not expect more than a 
small fraction of them to explode in the denser phases of the ISM. 
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7.3. Propeller-based Models for SGRs & AXPs 



The rapid spin-down rates, ages, clustered spin periods, and x-ray luminosities of 
AXPs and/or SGRs can all be explained by models involving the propeller effect as the 
dominant spin-down torque ( [van Paradijs et al. 1995| ; Alpar 2000; Chatterjee, Hernquist & 
Narayan 2000; Chatterjee Sz Hernquist 2000). In the context of these models, SGR bursts 
can be explained in terms of solid body accretion flHartwitt fc Salpeter 1973| ; [Tremaine fc 



Zytkow 1986|) or crustal instabilities ( Blaes et al. 1990]) . Because there are upper limits on 



the masses of the possible binary companions of at least some of the AXPs (e.g. |Mereghe"ttT 



Israel, fc Stella 1998| ), the accreted material is probably ejecta from the neutron star's own 



supernova explosion (|van Paradijs et al. 1995 ; |Corbet et al. 1995|) . In addition, since some 



of the SGRs and AXPs are located outside their apparent supernova remnants, the accreted 
material must form an accretion disk which can store the angular momentum of the accreted 
material and spin-down the neutron star on time scales of 1 — 10 kyr. Upper limits on the 
disk emission from the nearby AXP 2259+58.6 ( |Goe fc Pightling 1998| |Hulleman et al. 



2000 ) may rule out a standard hydrogen accretion disk around this AXP, but a disk may 
still exist given the uncertainties in the spectra of disks formed from metal-rich supernova 
ejecta and dust (e.g. [Perna, Hernquist, fc Narayan 2000|) . In addition, the noisy spin-down 



of SGR 1806-20 flWoods et al. 2000|) and the 6.4 keV emission feature flStrohmayer fc 



Ibrahim 2000|) observed from an SGR 1900+14 burst are both consistent with accretion 



disks around these objects. There are two possible scenarios involving accretion from the 
ejecta of the supernova explosion which produces the neutron star: fallback disk accretion, 
"pushback" disk accretion, and accretion involving high velocity neutron stars. 

It was recently proposed that AXPs may be formed from neutron stars accreting 
material from "fallback" accretion disks (Chatterjee, Hernquist & Narayan 2000; Chatterjee 
& Hernquist 2000). These disks may be formed from ~ 0.001 - 0.1M o QMichel 1988) ; |Lin] 
|Woosley, fc Bodenheimer 1991| ) of inner ejecta material < 2 hours after the initial core 
collapse in a type II supernova explosion (e.g. |Chevalier 1989| ; |Woosley fc Weaver 1995| ). 



Since a total accreted mass of only 10 M Q is required to explain the spin-down of SGRs 
and AXPs via the propeller mechanism (e.g. Alpar 2000; Chatterjee, Hernquist & Narayan 
2000; Chatterjee & Hernquist 2000), only a very small fraction (~ 10 -5 ) of the fallback 
material must be accreted into a disk for this model to explain spin-down rates and ages of 
the SGRS and AXPs. Formation of such a fallback disk is limited to < 7 days after the core 
collapse because of heating of the ejecta by 56 Ni ( |Chevalier 1989| ). This is long before the 
remnant feels the external environment early in the Sedov phase, and therefore formation 
of an early fallback disk may not be compatible with the evidence for dense SGR and AXP 
progenitor environments. 

A model involving fallback disks which form later in the evolution of the SNR, 
however, could be consistent with the evidence of dense SGR and AXP environments. In 
particular, the expansion of slow-moving ejecta can be reversed by the Sedov phase reverse 
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shock resulting from the interaction between the blastwave and the external medium (e.g. 
McKee 1974j ; [Truelove fc McKee 1999|) , and the subsequent implosion could result in the 



formation of a "pushback" disk - so named because the SNR ejecta is pushed-back onto 
the star because of the interaction with the dense environment surrounding the progenitor. 
Pushback disk formation would probably be most likely for: 1) neutron star progenitors 
which experienced a large amount of mass loss prior to supernova, and 2) progenitors 
in the dense ISM, which would provide the necessary pressure to confine the wind mass 
near the star. This is exemplified by the circumstellar environments surrounding the 
progenitor of SN 1987A, which is surrounded by n ~ 10 2 — 10 3 in wind material and HII 
gas (e.g. |Chevalier fc Dwarkadas 1995|) . The pushback process begins at the "reversal time' : 



t rev ~ 0.5 — 1.0 kyr ( Truelove fc McKee 1999| ), and the pushback mass is ~ O.4M given a 
total ejecta mass of 10M Q . As with the fallback disks discussed above, the formation of a 
pushback disk from only a small fraction of this matter would be required to explain the 
spin-down of SGRs and AXPs. Since this later fallback occurs after the majority of 56 Ni 
and 56 Co decays, the formation of the accretion disk around the neutron star would not be 
limited by radioactive heating. Such a model may be a plausible explanation for the dense 
SGR and AXP environments, and needs to be explored in more detail. 

The final propeller-based scenario for SGRs and AXPs involves high velocity neutron 
stars (HVNSs) capturing disk material from co-moving supernova ejecta, as first suggested 
Van Paradijs, Taam, and van den Heuvel (1995). Although the exact mechanism by which 
neutron stars are given substantial "kick" velocities at birth is not known, observations 
show that the kick velocities exceed 500 km s _1 in approximately 20% of all neutron stars 
( Cordes fc Chernoff 1998| ). In addition, this kick velocity appears to be independent of the 



dipole moment of the neutron star, as indicated by population studies of isolated radio 
pulsars ( Cordes fc Chernoff 1998 ; Pcshpandc, Ramachandran, fc Radhakrishnan 1999 ) and 



observations of extremely high velocity stars with canonical neutron star magnetic fields 
of ~ 10 12 G (e.g. PSR B2224+65: |Romani, Cordes, fc Yadigaroglu 19971 ). As implied 
by the ratio 6*/6snr listed in Table 1, many of the SGR/AXP positions are significantly 
displaced from the apparent centers of their associated SNRs. These displacements imply 
that the SGR/AXPs may have systematically large transverse velocities, although there 
is considerable uncertainty in the actual velocities, mainly because of uncertainties in the 
ages of the associated remnants. In addition, the actual space velocities of the SGR/AXPs 
are larger by an unknown factor dependent on the viewing angle. It has been estimated 
( |van Paradijs et al. 1995| ) that a 1O _4 M accretion disk may be acquired by a high velocity 



neutron star as it moves through nearly co-moving supernova ejecta. However detailed 
calculations of the time-dependent accretion rate for the range of ISM densities, progenitor 
mass loss parameters, and neutron star magnetic fields, initial spin periods, and velocities 
are needed to properly constrain such a model. This is beyond the scope of the present 
paper, and will be left for future work. 



Additional evidence in favor of the pushback disk and HVNS models for SGRs and 
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AXPs may be provided by the observed number of these sources. For the pushback disk 
model, the number of SGRs and AXPs less than t years old is given by iV = r b f w f env t, where 
rf, is the Galactic neutron star birthrate, f env is the fraction of neutron star progenitors in 
the warm dense ISM, and f w is the fraction of neutron star progenitors which experience 
mass loss sufficient to form a pushback disk in dense ISM environments. Since the rate of 
mass loss from stellar winds is an increasing function of the initial main sequence mass, the 
fraction f w can be estimated by considering the minimum stellar mass which undergoes 
pre-supernova mass loss sufficient to decelerate its supernova ejecta at early times. To first 
order, this should occur when the total wind mass emitted during the progenitor's life 
is approximately equal to the supernova ejecta mass. From the solar metallicity stellar 
models of Schaller et al. (1992) - which include mass loss - this occurs for stars of initial 
main sequence masses greater than M min ~ 27 M Q . Using M min and a Salpeter IMF with 
a maximum and minimum neutron star progenitor mass of 4OM and 8M , respectively 
( |Woosley fc Weaver 1995| ), we obtain f w ~ 0.1. Assuming f env < 0.2 (as discussed in §2) 
and r b ~ 1/40 yr" 1 ( |van den Bergh fc McClure 1994T ), we obtain N < 15 SGRs and AXPs 
with ages less than t = 30 kyr. This is similar to the observed number (12) of these sources, 
which can be taken as evidence supporting the pushback disk model. A similar calculation 
is possible for the HVNS model. In this case, the expected number of SGRs and AXPs less 
than t years old is given by iV = rf,fh v f e nvt, where fhv of neutron stars with high space 
velocities and r& and f env are defined as before. Using the same values of f env and r&, and 
assuming fh v ~ 0.2 (velocities > 500 km s -1 ; |Uordes fc Chernoff 1998| ), yields N ~ 30 
expected SGR and AXP sources with ages less than 30 kyr for the HVNS model. Again, 
these numbers are in the right ballpark for the observed numbers of SGRs and AXPs. 



8. Summary 



We have shown that soft gamma-ray repeaters (SGRs) and anomalous x-ray pulsars 
(AXPs) are born in regions of the interstellar medium which are denser than the 
environments typical of young neutron stars. This suggests that the development of SGRs 
and AXPs may be related to their environments, and we examine the implications of this 
on magnetar and propeller-based models for SGRs and AXPs. The evidence of dense 
progenitor environments would be consistent with the magnetar model only if magnetars are 
born exclusively in dense environments, which does not appear to be the case if magnetars 
form only from the most massive stars. Propeller-based models for SGRs and AXPs 
involving the formation of accretion disks from supernova ejecta appear to be consistent 
with the evidence for dense progenitor environments since these environments may induce 
the formation of such disks. This may occur in two ways. Pushback disks may be formed 
from the infall of the innermost ejecta, pushed back towards the neutron stars by prompt 
reverse shocks from the interactions of the expanding remnants with massive progenitor 
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winds confined close to the stars by dense surrounding gas - producing rapid deceleration 
of the expanding ejecta and strong prompt reverse shocks (Truelove & McKee 1999). Fossil 
disks may also form around high velocity neutron stars accreting from nearly co-moving 
supernova ejecta, slowed by the strong prompt reverse shocks in such dense environments 
(van Paradijs et al. 1995). 
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A. Notes On Individual SGRs & AXPs 



Here we discuss and reference the SNR distances, radii, and ages listed in Tables 1 
and 2. The SNRs associated with all but three of the SGRs and AXPs have kinematic 
distances derived from associations with objects (e.g. HII regions or molecular clouds) 
having known distances. The distances of the other three SNRs - G29. 6+0.1, G346.6-0.2, 
and G42.8+0.6 - are derived from consideration of Galactic spiral arms along the line of 
sight. The remnants Kes 73, N49, W28, and CTB 109 have published age estimates based 
on measurements of shock velocities and/or x-ray temperatures, and for these remnants 
we adopt an age range corresponding to the minimum and maximum values quoted in the 
literature. The rest of the associated remnants have no published age estimates, other 
than those determined by assuming Sedov expansion and distances determined from £ — D 
relations, which are biased toward denser phases of the ISM. For these SNRs we adopt a 
broad range of possible ages. We take a lower age limit of U ow = max[tf e ,t v ], where tf e and 
t v are the minimum ages derived by assuming free expansion of the SNR (at 10 4 km s _1 ) 
and a maximum SGR/AXP transverse velocity of 2000 km s _1 , respectively. For an upper 
age limit, we take 30 kyr, which is the maximum estimated age for a SNR/radio pulsar 
association in Table 3. Finally, the hydrogen column densities (iV#, listed in Table 2) are 
determined from the best-fit spectral models for the x-ray spectra of the SGRs/AXPs. 

AXP 1841-045: This AXP is coincident with the SNR G27.4+0.0 (Kes 73), which 
is estimated from the x-ray temperature to be 0.5 — 2.5 kyr old ( [Helfand et al. 1991 



Gotthelf fc Vasisht 1997|) . The diameter of the radio remnant is ~ 4' ([Kriss et al. 1985|) 



and the distance from HI absorption is estimated to be 6.0 — 7.5 kpc flSanbonmatsu fc 



Helfand 1992p . A power law fit to the x-ray spectrum of the AXP yields a best-fit value of 



Njj = (2.7 - 3.4) x 10 22 cm" 2 QGotthelf fc Vasisht 1997) ), which is consistent with such a 
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distance (see Figure 2), and with the Njj of 
1997 ) of the x-ray spectrum of the SNR. 



2 x 10 22 cm- 2 found (e.g. |Vasisht fc Gottheli 



SGR 0526-66: The error box of the SGR lies within the supernova remnant N49 ([Evans 



et al. 19"80| ; |Cline et al. 1980|) , located in the LMC at a distance of 49 — 55 kpc ([Feast 



1991| ; |Capccioli et al. 1990|) and interacting with a molecular cloud ([Hughes, Bronfman, & 



Nyman 1989j ). The diameter of N49 is ~ 1', and its estimated age is 5 — 16 kyr ( |Vancura et 



al. 1992| ; |Shull 1983| ) 



AXP 2259+58.6: The AXP lies within the ~ 30' diameter ( [Hughes et al. 1984Q SNR 
CTB 109 (G109. 2-1.0). The distance to the remnant, as determined by spectroscopy of 
stars in nearby HII regions, is in the range 3.6 — 5.5 kpc ([Rho fc Petre 1997|) . The estimated 
remnant age is 3.0 — 17.0 kyr ( |Parmar et al. 1998[ [Hughes, Harten, fc van den Burgh 
1981|) , and there is evidence for interaction between CTB 109 and surrounding molecular 



clouds (e.g. I'latematsu et al. 1987| ). The best-fit model to the x-ray spectrum of the AXP 



yields N H - 
N H = (0.9 



(0.9 — 1.0) x 10 22 cm 2 flParmar et al. 1998| ), which agrees with the range of 



1.1) x 10 22 cm 2 found for the x-ray spectra from various parts of the remnant 



( |Rho fc Petre 1997Q . 



AXP 1845-0258: This AXP is lies within the 5' diameter remnant G29.6+0.1 ( Gaensler 
et al. 1999| ). There are no kinematic distances, and the distance estimated from Galactic 



structure (see Figure 1) is 9.0 — 13.5 kpc. This is consistent (see Figure 2) with the AXP 

(9 ± 1) X 10 22 cm" 2 (|l'orii et al. 1998|). There are no reliable 



x-ray absorption of iV# = ± ij x iir~ cm 
SNR age estimates, so we adopt an SNR age range of 0.6 — 30.0 kyr by taking the limits 
of free expansion and the maximum detected SNR age from the radio pulsar sample. The 
chance association probability between the AXP and the SNR is estimated to be 1.6 x 10~ 3 
( IGaensler et al. 1999|) . 

SGR 1627-41: This SGR was first associated with the SNR G337.0-0.1 by Hurley et 
al. (1999c). The distance to this remnant is estimated to be 11.0 ± 0.3 kpc flCorbel et al. 



19991) , based on an association with a giant molecular cloud within the star-forming region 
CTB 33. The angular diameter of SNR G337.0-0.1 is ~ 1.5' QCorbel et al. 1999Q . The SNR 
age is unknown, so we adopt an age of 2.6 — 30.0 kyr. The probability that the SGR/SNR 
association is spurious was estimated to be ~ 5% ( |Smith, Bradt, fc Levine 1999| ). A 
marginal detection of a 6.4 s pulsation from the SGR was reported ( [Woods et al. 1999b|) 
but not confirmed by subsequent observations ( [Hurley et al. 2000b[ ). The x-ray spectrum of 
the SGR is equally well-fit by power law, blackbody, and thermal bremsstrahlung functions 
( |Hurley et al. 2000 b[ ), and the mean Nh for the power law spectral fit is (7.4 ± 0.6) x 10 22 



cm 2 ( |Woods et al. 1999b ; Hurley et al. 200"0b" ), which is consistent with the SNR distance 



(see Figure 2). 

SGR 1801-23: The most recently discovered SGR, this source has a long and thin 
error box (|Cline et al. 2000|) which passes through the center of SNR W28 (G6.4-0.1). W28 
is associated with OH masers in a molecular cloud whose distance rules out its previously 
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suggested association with PSR B1758-23 ( |Claussen et al. 1997|) . W28 has a radio diameter 
of ~ 40' (|Andrews et al. 1983Q , an estimated age of 2.4 — 30.0 kyr QLong et al. 199 1| ), and 



is located at a distance of 1.2 — 3.0 kpc flClark fc Caswell 1976| ; poudis 1976| ). The low 



energy absorption in the x-ray spectrum of W28 is Nh = (0.5 ± 0.1) x 10 22 cm 2 (|Rho fc 



Petre 1998|) , but there is no Nh estimate for the SGR since its burst spectra were were not 



measured to low enough energy to determine a value, and the persistent low energy x-ray 
counterpart of the source has not yet been identified. 



AXP 1709-40: This AXP lies ~ 8.5' form the center of the ~ 10' diameter flWhiteoak 



fe Green 1996| ) remnant G346.6-0.2. The chance association probability between the AXP 



and this remnant is ~ 0.1, from the method of §3. There are no published age estimates 
for the SNR, so we derive and age range of 3.6 — 30.0 kyr from the transverse velocity 
and maximum SNR age limits. There is no kinematic distance estimate to the SNR, so 
we adopt a distance of 3 — 5 kpc, which agrees with Galactic structure arguments placing 
it in or near spiral arms. This distance is also consistent (Figure 2) with the AXP x-ray 
absorption of N H = (1.81 ± 0.07) x 10 22 cm" 2 ( [Sugizaki et al. 1997|) . 



SGR 1806-20: This SGR lies within the ~ 7.5' diameter ( |Kulkarni et al. 1994|) 



remnant G10.0-0.3, which is located at a distance of 14.5 ± 1.4 kpc as determined from CO 
line observations of maser-associated molecular clouds flCorbel et al. 1997 ). This distance 



is consistent (Figure 2) with the large column density Nh = (6.0 ± 0.2) x 10 22 cm -2 found 
( (Sonobe et al. 19"94| ) for the x-ray spectrum of the SGR. The chance SGR/SNR association 
probability has been estimated at ~ 5 x 10~ 3 ([Kulkarni fc Frail 1993|) . The age is unknown, 
so we adopt an age of 3.5 — 30.0 kyr. This SGR and SNR may be associated with a compact 
cluster of stars ( [Fuchs et al. 1999"!) . 



SGR 1900+14: This SGR has been associated (e.g. [Hurley et al. 1999a| ) with the 24' 



diameter remnant SNR G42.8+0.6 ( [Fiirst et al. 1987|) from whose center it is displaced by 
17'. No kinematic distances to the SNR are available, so we adopt a distance of 3 — 9 kpc 
from association with a Galactic spiral arm (Figure 1). This distance is also consistent (see 
Figure 2) with the mean value of the SGR x-ray absorption for the best-fit spectral model: 
N H = (2.1 ± 0.2) x 10 22 cm~ 2 flWoods et al. 1999~a| ) . The age of the remnant is unknown, 



so we adopt an age range of 9.6 — 30.0 kyr. The chance probability of the SGR/SNR 
association by the method of §3 is ~ 0.1. This SGR has also been associated ( |Vrba et al.| 
|2000| ) with a compact star cluster at a distant of anywhere from 5 to 15 kpc, although the 
relatively low value of Nh from the absorption in the SGR x-ray spectrum (|Hurley et al 
|1999b| ; |Woods et al. 1999a| ) favors the near end of the range. 



AXP 1048-5937: This AXP is situated near the ~ 25' diameter remnant G287.8-0.5 
flJones 1973j ) associated with the Carina star-forming region at a distance of 2.5 — 2.8 kpc 
( [The fc Vleeming 1971| ; |Seward fc Chlebowski 1982|) . The geometrical chance association 



probability (§3) between the AXP and the remnant is ~ 0.16. There are no age estimates of 
the SNR and therefore we derive an age of 9.8 — 30.0 kyr from the transverse velocity limit 
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and the maximum detectable SNR age. The best-fit spectral model of the AXP emission 
yields the low value of N H = (0.45 ± 0.10) x 10 22 cm -2 QOosterbroeck et al. 1998Q , which is 
consistent (e.g. Figure 2) with the SNR/Carina distance. 



AXP 0720-3125: Despite its low x-ray luminosity (~5x 10 31 ergs s 1 , |Haberl et al 



1997|) , this 8.39 s pulsar has been included in the AXP sample due to its spin period, lack 
of an optical counterpart ( [Motch fc Haberl 1998| ), and fast spin-down rate (~ 3 x 10~ ] i 



Haberl et al. 1 9 97|) . Although there is no adjacent SNR in the catalog of Green (2000), the 



low x-ray absorption indicates a distance of only 0.10 ± 0.02 kpc ( [Haberl et al. 1997| ) to the 
source, which makes the detection of its associated SNR unlikely (see §3) if its age is as 
great (~ 40 kyr) as that indicated by its spin down rate. 



AXP 0142+615: This apparently old (possibly ~ 60 kyr) and well-studied (e.g. |White 



et al. 1996j ; [Israel, Mcrcghctti, fc Stella 1994| and references therein) AXP has no remnant 



within ~ 1° of it in the Green (2000) catalog. But there is evidence ( |White et al. 199~6|) 



that this source is in or behind a giant molecular cloud. Given the rapid evolution (e.g. 
l'ruelove fc McKee 1999| ) and possible self-absorption ( [Reynolds 1988| ) of the radio emission 



of SNRs in dense environments, the associated SNR may have faded below detectability, 
and therefore we are unable to constrain the radius of the associated remnant. 
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Fig. 1.— The positions of the soft gamma-ray repeaters (SGR) and anomalous x-ray 
pulsars (AXP) and their associated supernova shell remnants (SNR) projected onto a model 
( Taylor fc Cordes 1993| ) of the Galactic Plane distribution of free electrons. The position 
of the Galactic Center and the Sun are marked by the plus ( "+" ) and circle ( "o" ) symbols, 
respectively, and the line of sight to each source is shown by a dashed line. The estimated 
distances to the sources, from Galactic kinematic arguments or association with molecular 
clouds or spiral arms, are shown by the heavy solid lines. 
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Fig. 2. — The inferred dispersion measure flTaylor fc Cordes 1993| ) towards each of the 
SNRs, given the range of assumed distances, plotted against the observed N H values for the 
x-ray absorption of their associated SGRs and AXPs. As shown by the dotted lines, the data 
are consistent with the typical range of Galactic electron to hydrogen ratios (e/H) ( |Spitzer 
1978 ). This supports the SGR/AXP-SNR associations by showing that the SNR distances 
are consistent with the SGR and AXP distances. 
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SNR Age (yr) 

Fig. 3. — The radius of the SGR and AXP supernova remnant shells as a function of their 
age. The solid lines denote SNR expansion trajectories in the free expansion, Sedov and 
radiative phases (separated by dashed lines), according to Shull, Fesen, & Saken (1989), for 
an assumed supernova ejecta energy of 10 51 ergs in a wide range of ISM densities. The dotted 
lines denote the tracks of neutron stars born at the origin of the supernova explosion with 
varying space velocities. The data show that these objects are unusual in that they are all 
preferentially formed in the denser (> 0.1 H cm -3 ) phases of the interstellar medium (ISM), 
where <20% of all neutron-forming supernovae occur, as determined from observations of 
OB associations and Galactic supernova remnants. 
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Fig. 4. — Same as Figure 3, except for young radio pulsars with timing ages less than 32 
kyr. The remnants for these pulsars reside primarily in the diffuse phase of the ISM, as 
expected from observations of OB associations and Galactic supernova remnants. 
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Table 1. The Supernova Remnants of SGRs and AXPs 



Object 


Period 


SNR 


Age a 


Rad. b 


0* 

SSNR 






(s) 




(kyr) 


(pc) 


AXP 


1841-045 


11.8 


Kes 73 


1.5 + 1.0 


3.5-4.4 


0.1 


SGR 


0526-66 


8.0 


N49 


10 + 5 


7.1-8.0 


0.8 


AXP 


2259+58.6 


6.98 


CTB 109 


10 + 7 


16. - 24. 


0.2 


AXP 


1845-0258 


6.97 


G29.6+0.1 


15 + 14 


6.5-9.8 


0.1 


SGR 


1627-41 


6.4? 


G337.0-0.1 


16 + 14 


2.3 - 2.5 


2.3 


SGR 


1801-23 




W28 c 


16 + 14 


7.0-17.5 


0.1 


AXP 


1709-40 


11.0 


G346.6-0.2 C 


17+13 


4.4 - 7.3 


1.7 


SGR 


1806-20 


7.47 


G10.0-0.3 


17+13 


14. - 17. 


0.5 


SGR 


1900+14 


5.16 


G42.8+0.6 


20 + 10 


11. - 31. 


1.4 


AXP 


1048-5937 


6.45 


G287.8-0.5 C 


20 + 10 


9.1 - 10.2 


2.2 


AXP 


0720-3125 


8.39 


_f 


40 d 






AXP 


0142+615 


8.69 


_e 


60 d 







a SNR age (see Appendix notes) 

b Radius of radio shell (see Appendix and distances in Table 2) 
c "Tentative" remnant identification (see text) 
d MDR timing age since there is no identified SNR 
e In/behind molecular cloud (no identified remnant) 
f Too close to identify radio remnant 
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Table 2. Distance Measures of Galactic SGR and AXP and Associated SNRs a 



SNR 


I 


b 


Dist. b 


DM C 


SGR/AXP 






(deg.) 


(deg.) 


(kpc) 


(pc cm -3 ) 




(10 22 


cm 2 ) 


Kes 73 


27.386 


-0.006 


6.0-7.5 


440. 


- 620. 


AXP 1841-045 


2.7 


-3.4 


CTB 109 


109.093 


-0.993 


3.6-5.5 


100. 


- 135. 


AXP 2259+58.6 


0.9 


- 1.0 


G29.6+0.1 


29.679 


-0.109 


9.0 - 13.5 d 


860 


- 1200 


AXP 1845-0258 


8.0 - 


- 10.0 


G337.0-0.1 


336.968 


-0.111 


10.7-11.3 


1020. 


- 1070. 


SGR 1627-41 


6.8 


-8.0 


G346.6-0.2 


346.482 


+0.035 


3.0 -5.0 d 


120. 


- 340. 


AXP 1709-40 


1.7 


- 1.9 


G10.0-0.3 


9.996 


-0.242 


13. - 16. 


1070. 


- 1200. 


SGR 1806-20 


5.8 


-6.2 


G42.8+0.6 


43.021 


+0.766 


3.0-9.0 d 


70 


-460. 


SGR 1900+14 


1.9 


-2.3 


G287.8-0.5 


288.257 


-0.518 


2.5-2.8 


90. 


- 120. 


AXP 1048-5937 


0.4 


-0.6 



a See Appendix and references therein 
b Kinematic distance unless otherwise noted 

c Dispersion measures calculated from the model of Taylor & Cordes (1993) 
d Adopted distance (see text) 
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Table 3. Young Radio Pulsars 



Pulsar 


Period 


Age a 


Dist. b 


SNR 


Rad. c 


Reference 




(s) 


(kyr) 


(kpc) 




(pc) 




T0534+2200 


0.033 


1.3 


1.5 — 2.5 




> 17 


1 






1 6 


3 ^ — ^ 3 


IV-LkJ-LJ. -LtJ OA 


_1_ ' J _ _ 


9 


J0540-6919 


0.050 


1.7 


49.0 -55.0 


SNR 0540-693 


8-9 


3 


J1614-5047 


0.232 


7.5 


3.7-11.0 




> 30 


4 


J1617-5055 


0.069 


8.1 


6.1-6.9 




> 30 


5 


J0835-4510 


0.089 


11 


0.4-0.6 


Vela XYZ 


15-23 


11 


J1341-6220 


0.193 


12 


4.0 - 13.0 


G308.8-0.1 


15-48 


7 


J1801-2451 


0.125 


39 - 170 


2.3-6.9 


G5.4-1.2 


10-30 


8 


J1803-2137 


0.134 


16 


2.9-4.9 




> 30 


8 


J1709-4428 


0.102 


17 


1.4-2.7 




> 30 


2,10 


J1856+0113 


0.267 


20 


2.7-3.9 


W44 


12 - 17 


9 


J1048-5832 


0.124 


20 


2.2-3.7 




> 30 




J1740-3015 


0.607 


21 


2.5-4.1 




> 30 


6 


J1826-1334 


0.101 


22 


3.1-5.2 




> 30 


6 


J1730-3350 


0.139 


26 


3.2-5.3 




> 30 




J1646-4346 


0.232 


32 


3.4-10.3 


G341. 2+0.9 


9-28 


10 



^Timing age (from [laylor, Manchester, fc Lyne 1993] ; |Taylor et al. 1995|) , except 



for J1801-2451 ( |Gaensler fc Frail 200q) 
b Distance from pulsar dispersion measure 
c Radius of SNR radio shell 

References. - (1) [Frail et al. T995| ; (2) Caswell et al. 1981; (3) Manchester et 
al. 1993; (4) |Johnston et al. 1995| ; (5) |Kaspi et al. T99l| (6) Braun et al. 1989; (7) 



Caswell et al. 1992| ; (8) |Frail, Kassim, fc Weiler 1994| ; (9) Giacani et al. 1997fc (10) 



brail, Goss, fc Whiteoak 1991 ; (H) |Green 1984 
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Table 4. Occurrences in the Warm and Hot Phases of the ISM 

Source WISM(%) a HISM(%) b 

ExtraGal SNe < 20 > 80 

Galactic SNR 10 ± 10 90 ± 10 

Young PSRs 31 ± 14 69 ±21 

AXP/SGRs 83 ±26 17 ±12 



a Percentage in the warm ISM (n > 0.01 cm 3 ) 
b Percentage in the hot ISM (n ~ 0.001 cm -3 ) 

Note. - ExtraGal SNe - from van Dyk et al. (1996) 
observations of 49 extragalactic Type II & Ib/c SNe, 
corrected for detection threshold as discussed in text; 
Galactic SNR - from Higdon & Lingenfelter (1980) analysis 
of Clark & Caswell (1976) catalog of Galactic SNRs; Young 
PSRs - youngest (< 30 kyr) pulsars, see Table 3 and Figure 
4; AXP/SGRs - see Table 1 and Figure 3. 



